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In 1929, Leo Szilard invented a feedback protocol[? ] in which a hypothetical intelligence called 
Maxwell's demon pumps heat from an isothermal environment and transduces it to work. After an 
intense controversy that lasted over eighty years; it was finally clarified that the demon's role does not 
contradict the second law of thermodynamics, implying that we can convert information to free energy 
in principle!? ????]. Nevertheless, experimental demonstration of this information-to-energy 
conversion has been elusive. Here, we demonstrate that a nonequilibrium feedback manipulation of a 
Brownian particle based on information about its location achieves a Szilard-type information-energy 
conversion. Under real-time feedback control, the particle climbs up a spiral-stairs-like potential ex- 
erted by an electric field and obtains free energy larger than the amount of work performed on it. 
This enables us to verify the generalized Jarzynski equality[? ], or a new fundamental principle of 
"information-heat engine" which converts information to energy by feedback control. 

To illustrate the basic idea of our feedback protocol, let us consider a microscopic particle on a spiral-stairs-like 
potential (Fig. \T}. We set the height of each step comparable to the thermal energy k B T, where fc B is the 
Boltzmann constant and T is the temperature. Subjected to thermal fluctuations, the particle jumps between steps 
stochastically. Although the particle sometimes jumps to an upper step, downward jumps along the gradient are 
more frequent than upward jumps. In this manner, on average, the particle falls down the stairs unless it is externally 
pushed up (Fig. \Vp). Now, let us consider the following feedback control. We measure the particle's position at 
regular intervals, and if an upward jump is observed, we place a block behind the particle to prevent subsequent 
downward jumps (Fig. QJ)). If this procedure is repeated, the particle is expected to climb up the stairs. Note that, 
in the ideal case, energy to place the block can be negligible; this implies that the particle can obtain free energy 
without any direct energy injection. In such a case, what drives the particle to climb up the stairs? This apparent 
contradiction to the second law of thermodynamics, epitomized by Maxwell's demon, inspired many physicists 
to generalize the principles of thermodynamics!? ? ? ]. It is now understood that the particle is driven by the 
"information" gained by the measurement of the particle's location[? ? ]. 

In microscopic systems, thermodynamic quantities such as work, heat, and internal energy do not remain constant 
but fluctuate!? ? ]■ l n fact stochastic violations of the second law have been observed[? ? ]; nonetheless, the second 
law still holds, on average, if the initial state is in the thermal equilibrium : (AF - W) < 0, where AF is the free 
energy difference between states, W the work performed on the system, and (•) the ensemble average. However, 
the feedback control enables us to selectively manipulate only those fluctuations that violate the second law such as 
upward jumps by using the information about the system[? ? ? ]. Our Gedankenexperiment shows that by employing 
feedback control, the information can be used as a resource for free energy. In fact, Szilard has developed a model 
that converts one bit of information about the system to fc B Tln2 of free energy or work [? ]. In other words, the 
second law is generalized!? ] as follows: 

{AF-W)<k B TL (1) 

Here I is the mutual information content obtained by measurements!? ? ] (see Methods). So far, the idea of a 
simple thermal rectification by feedback control has found applications such as the reduction of thermal noise[? ] 
and the rectification of an atomic current at low temperature[? ]. On the other hand, the Szilard-type Maxwell's 
demon enables us to evaluate both the input (utilized information content) and the output (obatined energy) of the 
feedback control and relate them operationally. Therefore, it has provided an ideal test-ground of information-energy 
conversion and played the crucial role in the foundation of thermodynamics. However, its experimental realization 
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has been elusive. In this experiment, we developed a new method to evaluate the information contents and 
thermodynamic quantities of feedback systems and demonstrate the Szilard-type information-energy conversion for 
the first time using the colloidal particle on a spiral-stair-like potential. 

A dimeric particle comprising polystyrene beads (diameter = 287 nm) was attached to the top glass surface of 
a chamber filled with a buffer solution(Fig. |2k). The particle was pinned at a single point by a linker molecule; it 
exhibited rotational Brownian motion. By using quadrant electrodes imprinted on the bottom glass plate, we imposed 
1 -MHz electric fields to simultaneously create periodic potentials and constant torque on the particle along the angle 
of rotation. By using this new method, a tilted periodic potential with an ideal sinusoidal shape for the particle can be 
realized, which is a realization of the spiral-stairs-like potential mentioned above (Fig. |£t>, see also SI). A feedback 
control was performed under a microscope by constructing a real-time feedback system including video capture, 
image analysis, potential modulation, and data storage. We repeated the following feedback cycle with a period of 
t = 44 ms and a minimum feedback delay of 1 .1 ms, as illustrated in Fig. [2fc. At t = 0, the particle's angular position 
is measured. If the particle is observed at the angular region indicated as "S," the potential is changed to the one 
with an opposite phase at t = e; otherwise, no action is taken. At t = r, the next cycle begins with the measurement 
of the angular position. The region S was chosen for its energy advantage; in the region S, the potential energy 
before switching is always higher than that after switching. In the case of small e, the particle is expected to be 
at rest around the region S just before the switching at t = e and then jump to the rightward well of the switched 
potential after the switching. On the other hand, for large e, the particle falls down in the well away from the region S 
before the switching. In this case, with a large probability, the particle jumps down to the leftward well of the switched 
potential after the switching. In this manner, the feedback delay s regulates the efficiency of the feedback control. 
Note that, since r = 44 ms is sufficiently larger than the relaxation time in each well (~10 ms) and smaller than the 
typical time to jump to neighbor wells (—1 s), each feedback cycle is supposed to be a transition between equilibrium 
states. 

In Fig. [3k, typical trajectories with the feedback control are shown. The trajectories are stepwise with a step size of 
90 degrees, which reflects the potential profile (Fig. |3p). We find that for small e, the particle rotates uni-directionally 
while climbing up the potential, whereas for large e, the particle goes down along the gradient. The rotation rate 
decreases monotonically with e, as expected (Fig. |3fc). 

We then focused on the energetics during a cycle. In Fig. [3jJ, we show the difference between the obtained 
free energy AF and the work performed on the particle by the switching, W, which is averaged over a cycle (see 
Methods). We find that (AF - W) > for small e\ this implies that the particle gains a net free energy larger than 
the work performed by absorbing heat beyond the conventional limitation of the second law of thermodynamics. For 
small e, the switching mostly occurs when the particle is in the region S. In such cases, the particle absorbs heat 
from an isothermal environment to reach the region S before the measurements at t = 0, then performs work to the 
electric field at the switching, and finally jumps to the rightward well after the switching (Fig. S7). Although such an 
event is not prohibited even if we randomly switch potentials without feedback control, it is typically an accidental 
and rare event in accordance with the second law of thermodynamics or the fluctuation theorem[? ? ? ]. However, 
the feedback control can increase the likelihood of occurrence of such an event. This is the crux of the control 
by Maxwell's demon. The resource of the excess free energy is the information obtained by the measurement. If 
the estimation error of the particle's angular position is negligible, the amount of information is characterized by 
the Shannon information content I [? ]. In this study, I = —phxp - (1 - p)ln(l - p), where p is the probability 
that the particle is observed in the region S (see Methods). As noted in (TJ, I can be converted to free energy 
of up to k B TI [? ]. In our system, for the shortest feedback delay (e = 1.1ms), p, I, and (AF - W) were 0.059, 
0.22, and 0.062 fc B T, respectively (Fig. S9). This gives the efficiency of the information-energy conversion as, 
(AF - W)/k B TI = 28%. The 100% efficiency can be achieved by quasi-static information heat engines such as the 
Szilard engine [? ]. 
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Although the second law concerns only the average, or the first-order cumulant, of the stochastic quantity AF - W, 
Jarzynski pointed out that the second law naturally emerges as the first-order cumulant expansion of the following 
equality that involves AF - W to all orders [? ? ]: ( e ( AF ~ w ^ kBT ) = 1. Recently, the Jarzynski equality, which 
assumes a prescribed control scheme, was generalized to systems with a feedback control as follows[? ]: 

{e (AF~W)/k BT) = 7j (2) 

where 7 is an experimentally measurable quantity and defined as the sum of the probabilities that the time-reversed 
trajectories are observed under time-reversed protocols for all possible protocols (see Methods). From its definition, 
< 7 < 2 in our system. While I concerns the information obtained by the measurements, 7 quantifies how 
efficiently we utilize the obtained information for the feedback control. If we control the system perfectly and 
deterministically, a time-reversed trajectory is always realized under the time-reversed protocol; 7 then takes the 
maximum value. We repeated time-reversed cycles with/without switchings to obtain 7 with a period of 220 ms, 
which is sufficiently long to ensure that the initial state of the cycle is relaxed to equilibrium. Figure |4]shows that the 
conventional Jarzynski equality is violated in the presence of the feedback control. For large e where (AF - W) < 0, 
the second law holds on average, but the Jarzynski equality is violated. On the other hand, the generalized Jarzynski 
equality |2) holds over a broad range of e, showing that |2) expresses the effect of feedback control to all orders. 7 
appears to converge to 1 in the limit of infinite e; here, the angular position at the switching becomes independent 
of that at the measurement, and the conventional Jarzynski equality recovers. The validity of © verifies a new 
fundamental principle of an "information-heat engine", which converts information to free energy, in the all-order 
terms. 

We have demonstrated the information-heat engine, which converts information to free energy, as the first realiza- 
tion of a Szilard-type Maxwell's demon. Note that, since the energy converted from information is compensated for 
by the demon's energy cost to manipulate information!? ? ? ], the second law of thermodynamics is not violated 
when the total system including both the particle and demon is considered. In our system, the demon comprises 
of macroscopic devices such as computers, the microscopic device gains energy at the expense of the energy con- 
sumption of a macroscopic device. In other words, by using information as the energy-transferring "medium", this 
information-energy conversion can be utilized to transport energy to nanomachines [? ? ] even if it is not possible to 
drive them directly (Fig. S1). It is the future challenge to realize a nanoscale information-processing device such as 
an artificial molecular motor[? ? ], in which both the demon and controlled system are microscopic. 

METHODS 
Free energy and work. 

Each potential well separated by peaks was defined as a state (Fig. S5). The free energy of the state fc was calu- 
clated as F k = -k B T\n [J dxe _c/(;r)/ ' CBT ], where the integration was performed in the angular region corresponding 
to the state k. Since the shapes of all the wells are almost the same, the free energy difference between states is 
nearly equal to the difference of the potential energy of their local minima. The work performed on the particle W 
was calculated as the potential-energy change associated with the switching; the potential energy after the switching 
minus that before the switching. In cycles without switching, W = 0. 

Information content. 

For an event k with a probability of occurrence p(k), the Shannon information content associated with this event 
is defined as -lnp(fe). This definition leads to well-defined properties that the information content should satisfy[? 
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]. The average Shannon information content becomes I = -£) fc p(fc)]np(fc). Measurements usually accompany 
errors, which reduce the amount of information that can be utilized. Although I denotes the amount of the information 
embedded in the system, the mutual information content, I', denotes the amount of information that is obtained by 
the measurement [? ?]:/' = J2k, m P( m \ k )P( k ) ln P< kk)' > ■ wnere p(m\k) is the conditional probability that the outcome 
of the measurement is the m-th event when the fc-th event occurs actually. If the measurement is free from error, 
p(m\k) = 5k, m (Sk,m = 1 if k = m, and otherwise 0). In such case, I' = I. In the present experiment, we distinguished 
two events; the particle is observed in the region S or not with negligible measurement errors. Then, the (average) 
Shannon information content per cycle becomes the so-called binary entropy function: I = -plnp-(l-p) ln(l -p), 
where p is the probability that the particle is observed in the region S. 

Feedback efficacy. 

The feedback efficacy 7 is defined as the sum of the probabilities that the time-reversed trajectories are observed 
under time-reversed protocols for all possible protocols (Fig. S7). In the forward feedback cycle, we measure the 
particle's angular position at t = and (i) switched or (ii) did not switch the potential at t = s depending on the angular 
position. Corresponding time-reversed trajectories are that the particle is observed in the region (i) S at t = r after the 
switching at t = t - e or (ii) outside S without switching. Let the occurrence probabilities of time-reversed trajectories 
under timer-reversed protocolos be p sw and p ns , respectively. Then, 7 is 7 = p sw +p ns - From its definition, if there are 
m states to be distinguished (m = 2 in our experiment: whether the particle is in the region S or not), < 7 < m. We 
repeated time reversed cycles with/without switchings to obtain 7 with a period of 220 ms, which is sufficiently long 
to ensure that the initial state of the cycle is relaxed to equilibrium. 
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FIG. 1 : Schematic illustration of the experiment, a, A microscopic particle on a spiral-stairs-like potential with a step height 
comparable to k B T. The particle stochastically jumps between steps due to thermal fluctuations. Since the downward jumps 
along the gradient are more frequent than the upward ones, the particle falls down the stairs, on average, b, Feedback control. 
When an upward jump is observed, a block is placed behind the particle to prevent downward jumps. By repeating this cycle, the 
particle is expected to climb up the stairs without direct energy injection. 




FIG. 2: Experimental setup[? ? ]. a, A dimeric particle composed of particles with a diameter of 287 nm (Seradyn) was non- 
specifically attached to the top glass surface via a streptavidin linker coated on the particle's surface. The particle was pinned at 
a single point and exhibited rotational Brownian motion. To impose a tilted periodic potential on the particle, an elliptically rotating 
electric field (blue and gray curves) was induced by applying 1 MHz sinusoidal voltages on the quadrant electrodes patterned on 
the bottom glass surface. The particle was observed on a microscope equipped with a high-speed camera at a period of 1 .1 ms 
with an exposure time of 0.3 ms. Not to scale. See SI for details, b, Typical potentials with opposite phases to be switched in 
the feedback control. Potentials were measured from transition probabilities (see SI). The particle experienced a tilted periodic 
potential with a period of 180 degrees. The direction of the long axis of the elliptically rotating electric field corresponds to the 
local minima of the potential. By changing the direction of its axis, we inverted the phase of the potential. The height and slope 
were 3.05± 0.03 k B T and 1.13 ± 0.06 k B T/360° (meaniS.E., 7 particles), respectively, c, Feedback control. At t = 0, the 
particle's angular position is measured. If the particle is observed in the angular region indicated by "S," we switch the potential at 
t — e by inverting the phase of the potential (right). Otherwise, we do nothing (left). At t = t, the next cycle starts. The location 
of the region S is altered by the switching. Potential's wells correspond to the steps of the spiral stairs in Fig. Q] The switching of 
potentials corresponds to the placement of the block. 
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FIG. 3: Trajectories, mean velocities, and excess free energy under feedback control, a, Typical trajectories for different values of 
the feedback delay e. b, Magnified plot of the region indicated by a rectangle in (a). The particle rotates with steps with a step 
size of 90 degrees reflecting the profile of the potential. Red triangles indicate the timings of the switchings, c, Variation of the 
rotation rate with feedback delay e. The rotation rate is defined as positive when the particle climbs up the potential. Data of 
seven particles are averaged. Error bars indicate standard deviations, d, AF : the free energy difference between the initial and 
final states of the cycle; W : the amount of work performed on the particle by the switching calculated as the potential-energy 
change associated with the switching (see Methods). In cycles without switching, W — 0. (•) denotes the mean per cycle. In the 
shaded region, we obtain the excess free energy beyond the conventional limitation of the second law of thermodynamics. Error 
bars indicate standard deviations. 
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A. SCHEMATIC OF EXPERIMENT AND THE CONSISTENCY WITH THE SECOND LAW OF THERMODYNAMICS. 



We demonstrated that free energy is obtained by a feedback control using the information about the system; 
information is converted to free energy, as the first realization of Szilard-type Maxwell's demon[1]. 

Since the obtained free energy or work is compensated for by the demon's energy cost to manipulate information, 
it does not violate the second law of thermodynamics when the total system including both the particle and demon 
is considered. The demon consists of macroscopic devices in practice such as computers in our system; the micro- 
scopic device gains energy at the expense of the energy consumption of a macroscopic device (Fig. S3). In other 
words, by using information as the energy-transferring "medium", this information-energy conversion can be utilized 
to transport energy to nanomachines even if it is not possible to drive them directly. 
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Figure S 1 : Schematic of the information-energy conversion with a macroscopic demon and a microscopic system (Szilard 
engine[1] as an example). The Szilard-engine can achieve the 100% conversion rate from information to energy; 1 bit (Ln2 
nat in the natural logarithm) of information is converted to free energy or work of fc B Tln2 in the microscopic system at the expense 
of an energy consumption of k B T\n2 in the macroscopic demon. 



B. PREPARATION OF CHAMBER. 

Brownian particles in a buffer solution were prepared for the experiment by the following procedure. The cover 
slip (Matsunami, Japan) that served as the ceiling of the observation chamber was washed three times with distilled 
water. The bottom glass coated with chromium electrodes (Matsunami, Japan) was soaked in 10 N KOH over night 
and then washed with distilled water. The top and bottom glasses were separated by a plastic film with a thickness of 
about 12 /iim (Saran wrap, Asahi Kasei, Japan) and silicone grease. 6 p\ of A buffer (50 mM MOPS-K, 50 mM KCI, 2 
mM MgCb, 2 mg/ml heat-shocked Bovine serum albumin (Sigma-Aldrich, USA), pH6.9) was applied to the chamber 
to suppress the binding of particles to the glass surface. Then, 10 f.i\ of streptavidin-coated polystyrene beads with a 
diameter of 287 nm (Seradyn, USA) in A buffer was applied to the chamber. Some particles bound to the top glass 
surface. About 50 fi\ of B buffer (5 mM MOPS-K, 1 mM MgCI 2 , 1 mM KP i: pH6.9) was applied to wash out unbound 
particles in the chamber. All the experiments were performed in B buffer at a room temperature. 
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C. MICROSCOPY AND IMAGE ANALYSIS. 



The incidentally-dimeric particle pinned at the top glass surface was observed on a phase-contrast microscope 
(BX-51 Wl, Olympus, Japan) equipped with a high speed camera (IPX-210L, Imperx, USA) at a period of 1.1 ms with 
an exposure time of 0.3 ms. Real-time feedback system was constructed on the PXI (PCI extensions for Instrumen- 
tation) system implementing a computing module (PXI-8108, National Instruments, USA), a video capturing board 
(PXI-1428, National Instruments), and a digital-to-analog conversion board for the control of the electric potential 
across the electrodes patterned on the bottom glass plate (PXI-6221 , National Instruments). The minimum delay in 
this feedback control is thus 1 .1 ms which is much shorter than the time for the particle to perform rotational diffusion 
up to 90 degrees. The angle of the dimeric probe particle was estimated using an algorithm based on the principal 
component analysis. 



To provide a periodic potential and adverse torque to the Brownian particle, we utilized the induced polarization 
effect of the dielectric material subjected to an AC external electric field[2-4]. Four chromium electrodes on the 
bottom glass surface were connected to four power amplifiers (HSA4101, NF, Japan) that amplify AC (Alternative 
Current) voltages generated by function generators (WF1974, NF, Japan). We applied 1-MHz AC voltages with 
some phase shifts on the electrodes (Fig. S0. Let the phases on the electrodes to be 0, 5, 1 80, and 180 + 5 degrees. 
Three different types of potentials are imposed on the particle according to the value of 5 (Fig. S©. 

• If <5 = ±90 degrees, a rotating electric field is applied at the center of the electrodes (Fig. S[3k), where a dipole 
moment rotating at the same rate with some time delay appears on dielectric particles. The length of the time 
delay is determined by the complex dielectric constants of the particle and the solution. Due to the time delay, 
the particle experiences a torque from the electric field [2-4]. The torque is constant; it does not depend on the 
particle's motions but depends only on the time delay. Its magnitude is proportional to the particles's volume 
and the square of the amplitude of the applied AC voltages. Also, it nonlinearly depends on the frequency of 
the applied AC voltages. According to the sign of 5, the direction of the torque is changed. 

• If 5 = or 180 degrees, the electric field does not rotate but oscillates along an axis (Fig. Sj3fc>). In this case, the 
particle feels a periodic potential with a period of 180 degrees along the direction of rotation. The oscillating 
direction of the electric field corresponds to the angle of the potential's local minima. 

• Otherwise, an elliptically rotating electric field is applied (Fig. $3fc); it corresponds to the superposition of the 
above two cases (6 = and 90 degrees). In such case, the particle feels a tilted periodic potential. The gradient 
of the potential depends on the value of 6. According to the sign of 5, the direction of the torque is changed. 



D. ELLIPTICALLY-ROTATING ELECTRIC FIELD. 
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Figure S 2: (a) Quadrant electrodes. Not to scale, (b) Sinusoidal voltages applied on the electrodes. 
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We set 6 to be about 5 degrees (6 was tuned for each particle) to apply the elliptically rotating electric field. For 
feedback control, the potential was switched by modulating the phases of applied voltages from 0, 6, 180, and 180+6 
degrees to 0, 180-6, 180, and -6 degrees, respectively for four electrodes. 
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Figure S3: Three different types of the potential imposed on the particle according to the phase shift 5. (a) Rotating electric field, 
(b) Oscillating electric field, (c) Elliptically rotating electric field. 



E. ESTIMATION OF POTENTIAL ENERGY. 

The potential profiles were estimated by a method similar to Ref. [5] (Fig. $4^). We switched two potentials with 
opposite phases periodically with a period of 220 ms to sample data around not only the local minima but also the 
peaks. We divided the whole video frames to two sets according to the potential at each frame. For each set of 
frames, we performed the following procedures to estimate the potential profile. 360 degrees are divided to equally- 
spaced bins with a bin width of 3 degrees. We counted the number of video frames nii where the particle is in the z'th 
bin. We also counted the number of transitions n^, or the number of two-successive-frame where the particle in the 
i th bin moves to the /th bin at the next frame. The transition probabilities = n^/nii satisfy the detailed balance 
condition: 

7/;. 

= e -AUH/feT (S1) 

where A!i;_» y - = Uj-Ui is the difference of the potential energies between z'th and / th bins. We searched the set {U*} 
which minimizes a score function as follows: 

£ 2 ({U,}) = £ [AIT M - AU?_J 2 , 

<<; 

where AU[^. = /c B r[lnz(;j^i - InH^j], using the Powell's method 
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Figure S4: (a) Potential estimation, (b) Estimation of the potential energy for a simulated trajectory. Solid lines indicate the 
theoretical curves. Markers correspond to estimated values (bin width = 3 degrees). The slope is 1 fc B T/360°, the activation energy 
(peak-to-peak amplitude) 3 k B T, the rotational friction coefficient 0.3 pN-nm-s/rad 2 , the framerate 1000 Hz, and temperature 300 
K. Number of frames is 1 ,000,000. The potentials were switched every 200 frames. 

The validity of this method was verified for numerically-simulated trajectories based on a Langevin dynamics. 
Figure SHfc> shows that this algorithm estimates the potential energy for the trajectory generated by a Langevin 
simulation to a good extent. 

A typical number of video frames used for the potential estimation was about 1 ,000,000. 

F. STATE AND FREE ENERGY. 

Figure S[5] is the schematic to illustrate the definition of the state and how to calculate the free energy. 




state k 



Figure S5: Definitions of the states and free energies. The well separated by peaks is regarded as a state. A free energy is 
assigned to each state. 
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G. SCHEMATIC OF ENERGETICS. 



Schematic to illustrate the energetics in a cycle with a switching is shown in Fig. S[6j 
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Figure S6: Schematic to illustrate the energetics in a cycle with a switching. For small e, switchings occur mostly when the 
particle is in the region S. In this case, the particle absorbs heat from an isothermal environment to reach the region S before 
measurements at t = 0, performs extractable work to the electric field at the switching, and finally jumps to the rightward well after 
the switching. On the other hand, for large e, most of the energy at t = dissipates to the heat bath before switching. Then, the 
switching performs external work on the particle with lifting up the potential, which then dissipates as it falls down in the well. 



H. FEEDBACK EFFICACY. 



Figure $7]is the schematic to illustrate the feedback efficacy. 

a Forward b 



t = [n 

Measurement 



t = e 



t = z 



[ Switch"*) 

IT 



Reverse 



t = T 

Measurement 

t = X —£ 



t = 



Psw \ / P 



( Switch") 

; i Pns ~Y~ 



Figure S7: Forward and reversed cycle. (A) Forward feedback cycle. "S" and "N" indicate that the particle is observed in the 
region S or not, respectively. (B) Reversed cycle. Feedback control is not performed. 
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I. EFFICIENCY. 

Figure S[8]shows (a) the probability that the particle is observed in the region S, (b) the (average) Shannon infor- 
mation content, and (c) the efficiency to convert information to energy. 
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Figure S8: Efficiency of information-energy conversion, (a) The probability p that the particle is observed in the region S. (b) 
Shannon information content J = -p\np - (1 - p)ln(l - p). (c) Efficiency of the information-energy conversion calculated as 
<AF - W)/k B TI. Error bars indicate standard deviations. 
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